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We first define a series of AW interaction models ranging from very simple to fully realistic. We 
then present Green's function Monte Carlo calculations of light nuclei to show how nuclear spectra 
evolve as the nuclear forces are made increasingly sophisticated. We find that the absence of stable 
five- and eight-body nuclei depends crucially on the spin, isospin, and tensor components of the 
nuclear force. 

PACS numbers: PACS numbers: 21.10.-k, 21.45.-|-v, 21.60.Ka 
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A key feature of nuclear structure, of great importance 
to the universe as we know it, is the absence of stable 
five- or eight-body nuclei. This simple fact is crucial to 
both primordial and stellar nucleosynthesis. It leads to a 
universe whose baryonic content is dominated by hydro- 
gen and '*He, with trace amounts of deuterium, '^He, and 
^Li. It also enables stars like our sun to burn steadily for 
billions of years, allowing time for the evolution of life 
intelligent enough to wonder about such issues. 

In this Letter we demonstrate that the binding ener- 
gies, excitation structure, and relative stability of light 
nuclei, including the opening of the A = 5 and 8 mass 
gaps, are crucially dependent on the complicated struc- 
ture of the nuclear force. We do this by calculating the 
energy spectra of light nuclei using a variety of nuclear 
force models ranging from very simple to fully realistic, 
and observing how features of the experimental spectrum 
evolve with the sophistication of the force. Wc find that 
the spin-isospin and tensor forces present in long-range 
one-pion-exchange (OPE) are vital, which in turn may 
allow us to make a closer connection between nuclear 
structure and the underlying features of QCD ||]. 

Modern nucleon-nucleon {NN) potentials, such as the 
Argonne vis §], CD Bonn |], Reid93, Nijm I, and 
Nijm II j|] , fit over 4300 elastic NN scattering data with a 
« 1. These potentials are very complicated, including 
spin, isospin, tensor, spin-orbit, quadratic momentum- 
dependent, and charge-dependent terms, with ^-^40 pa- 
rameters adjusted to fit the data. Despite this sophis- 
tication, these potentials cannot reproduce the binding 
energy of few-body nuclei like and "^He without the 
assistance of a three-nucleon potential . Three- nucleon 
(NNN) potentials, such as the Tucson-Melbourne Ur- 
bana and Illinois |^ models, are also fairly compli- 
cated, depending on the positions, spins, and isospins 
of all three nucleons simultaneously. A combination of 
AW and NNN potentials, such as the Argonne vis and 
Illinois 2 (AV18/IL2), evaluated with exact Green's func- 
tion Monte Carlo (GFMC) many-body calculations, can 
describe the spectra of light nuclei very well |l^ . 

The AV18 potential contains a complete electromag- 
netic (EM) interaction and a strong interaction part 



which is a combination of OPE and remaining shorter- 
range phenomenology. The strong interaction part is 
written as a sum of 18 operator terms: 



p=l,18 



The first eight operators. 
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are the most important for fitting S- and P-wave AW 
data. The additional terms include six operators that are 
quadratic in L, three charge-dependent (CD) terms, and 
one charge-symmetry-breaking (CSB) term. The radial 
functions Vp{r) have parameters adjusted to fit the elastic 
pp and np scattering data for Eiab < 350 MeV, the nn 
scattering length, and the deuteron energy. 

The Illinois NNN potentials include a complete two- 
pion-exchange piece, three-pion rings, and a shorter- 
range phenomenological term: 
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The five Illinois models (IL1-IL5) each have ~3 parame- 
ters adjusted to fit the energies of 17 narrow states in 
A < 8 nuclei as evaluated in GFMC calculations [1). 
Subsequent calculations of an additional 10 states in 
A = 9,10 nuclei show that, without readjustment, the 
AV18/IL2 combination is able to reproduce 27 narrow 
states with an rms deviation of only 600 keV |l0[| . 

The AV18/IL2 Hamiltonian is the standard of com- 
parison for this Letter. We also present previously re- 
ported results for AV18 alone, and for a simplified, but 
still fairly reahstic, potential called AV8' ||ll|. The "8" 
designates the number of operator components, which in 
this case means those of Eq.(||). The standard Coulomb 
interaction between protons, Vci{pp), is retained, but 
smaller EM terms are omitted. The prime denotes that 
this potential is not a simple truncation of AV18, but 
a reprojection, which preserves the isoscalar average of 
the strong interaction in all S and P partial waves as 
well as in the '^Di wave and its coupling to ^Si. Conse- 
quently, the deuteron bound state is virtually identical to 
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that of AV18, except that the omission of the small EM 
terms alters the binding energy from the experimental 
value of 2.22 MeV to 2.24 MeV. Details of this reprojec- 
tion are given in Ref. Recently, the AV8' (without 
Coulomb) was used in a benchmark calculation of '^He by 
seven different methods, including GFMC, with excellent 
agreement between the various results | p^ . 

Here we define five new potentials, which are progres- 
sively simpler reprojections of AV8', designated AV6', 
AV4', AVX', AV2', and AVI' The reprojections 

preserve as many features of elastic NN scattering and 
the deuteron as feasible at each level. GFMC calcula- 
tions of A < 10 nuclei for these simpler models show how 
different features of the spectra correlate with specific 
elements of the forces. 

The AV6' is obtained by deleting the spin-orbit terms 
from AV8' and adjusting the potential to preserve the 
deuteron binding. The spin-orbit terms do not contribute 
to S-wave NN scattering, and are the smallest contrib- 
utors to the energy of ^He They also do not con- 
tribute to scattering in the ^Pi channel, but are impor- 
tant in differentiating between the '^Po.1,2 channels. Thus 
this model does not give a very good fit to NN scatter- 
ing data. To fix the deuteron, we choose to subtract a 
fraction of the AV8' spin-orbit potential's radial function 
from the central potential in the ST = 10 channel, ad- 
justing the coefficient to get an energy of 2.24 MeV: 

v-io{AV6') = v1q{AV8') - 0.3 v[l{AV8') . (4) 

This choice preserves the OPE potential, while the 
deuteron D-state and quadrupole moment are barely 
changed. Spin and isospin terms are projected from the 
Vgj, as in Eq. (30) of Ref. |^ , while tensor terms remain 
the same as in AV8'. 

The AV4' potential eliminates the tensor terms. The 
^So and ^Pi potentials are unaffected, but the coupling 
between ^Si and "^Di channels is gone and the '^Po,i,2 
channels deteriorate further. The central ST = 10 po- 
tential is again adjusted to fix the deuteron binding: 

vIoiAVA') = v1o{AV6') + 0.8735 vIq{AV6') , (5) 

but now there is no D-state and no quadrupole moment. 

Although many aspects of NN scattering have been 
sacrificed at the AV4' level, such a potential still allows 
us to differentiate between the four possible ST channels. 
Any further reduction in the operator structure sacrifices 
this feature. We consider three such simplifications: 1) 
AVX', where the operators are 1 and the space exchange 
operator P^j = -^{1 + ai ■ aj + n ■ r, -I- ■ a^n • r,); 2) 
AV2', with operators 1 and • r,; and 3) AVI', which is 
just a pure central force. AVX' allows one to differenti- 
ate between the spin-isospin weighted averages of S-wave 
and P-wave forces by setting v'^ + — ^(woi + wfo) and 

_ _ _L.^g^c^ _|_ y^^y In this case, the average of 
the S-waves gives a deuteron that is bound by only 0.43 



MeV, but the intrinsic repulsion in odd partial waves is 
retained. AV2' allows one to differentiate between ^S and 
^S potentials, analogous to the Malfiiet-Tjon (MT) I-III 
interaction |l^] , with the combinations v'^ + v'^ = v^i and 
v'^ — 3w'^ = ViQ. Finally, AVI' is just the average of ^S 
and '^S potentials, analogous to the MT V interaction. 
The AV2' preserves the correct deuteron binding of 2.24 
MeV, but the AVI' again has a deuteron bound by only 
0.43 MeV. While the MT interactions were intended only 
for use in s-shell {A < 4) nuclei, they have been used in 
larger systems acting either in all partial waves, or only 
in even partial waves Here we treat AVI', AV2', and 
AVX' as operators acting in all partial waves. 

In the independent-particle (IP) model, nuclear states 
can be characterized by quantum numbers '^^^^L[n] as 
well as , where L and S are the orbital and spin angular 
momenta, [n] designates the Young tableau of the spatial 
symmetry, and are the total angular momentum and 
parity |jl^ . Realistic interactions down to the AV6' level 
mix states of different ^^'^^L[n], but for the AV4' and 
simpler interactions there is no mixing of different values 
of S and essentially no mixing of L or [n] . 

Our many-body calculations are made with the GFMC 
method, details of which may be found in Refs. 1,0 |l|: 
. The GFMC method is in principle exact for the AV8' 
and simpler potentials, while for the AVIS and AV18/IL2 
models some small parts of the interaction have to be 
calculated perturbatively. We believe the calculation of 
binding energies for AV18 and AV18/IL2 is accurate to 
1-2%, and better for the simpler models. Results for 26 
'^^~^^L[n] states in nuclei ranging from ^He to "'^"B are 
shown in Fig. |l]. The Hamiltonians are, from left to right 
for each nucleus, AVI', AV2', AVX', AV4', and AV6'. 
Results for 25 J'^ states are shown in Fig. ^ for AV6', 
AV8', AV18, AV18/IL2, and experiment @. 

The simple AVI' and AV2' interactions approximately 
reproduce the energies of s-shell nuclei ('^H is overbound 
by ^ 0.6 MeV). However, every additional nucleon sig- 
nificantly increases the binding, so p-shell nuclei are pro- 
gressively more bound and no mass gaps appear. With 
such forces, nuclear matter will not saturate until the re- 
pulsive cores of the interaction start to overlap at many 
times the empirical density. For AVI', the strong in- 
teraction does not differentiate between different isospin 
states, so ^He and ^Li would have the same energy except 
that the Coulomb interaction makes ^Li less bound. Con- 
sequently, the /^-stable nuclei in the 4 < A < 10 regime 
would all be isotopes of helium. The AV2' avoids this 
particular problem by preserving separately the deuteron 
binding and ^S scattering; it also improves the saturation 
behavior slightly. Both models have the curious feature 
that the spectrum is reversed in order, i.e., the lowest p- 
shell eigenstates are the ones that are spatially least sym- 
metric. This lowers the energy by reducing the amount of 
overlap of the repulsive potential cores in the wave func- 
tion. Thus the ground state for ^Be is a '^P[211] state. 
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FIG. 1: Nuclear energy levels for the simpler potential models; dashed lines show breakup thresholds. 



which has degenerate spin of J"^ = 0+, 1+, 2+. The spec- 
tra are also very compressed compared to experiment. 

The AVX' and AV4' overcome many of the limitations 
of the simplest models by preserving the difference be- 
tween attractive even- and repulsive odd-partial waves. 
Both provide significant saturation, particularly the fea- 
ture that A — 5 nuclei are unstable. However, the A — 8 
mass gap is a more subtle effect, since both these models 
predict ^Be to have slightly more than twice the binding 
of ''He. The lowest states are the spatially most symmet- 
ric, so *Be now has a ^S[4] = 0+ ground state. The 
spectrum is also less compressed. Because the AVX' does 
not differentiate between 'S and '^S interactions, it shares 
the failing of AVI' in having ^He more bound than ^Li, 
but due to the correct ordering of spatial symmetries, 
®Be is much more bound than ^He. 

The tensor forces in AV6' provide significant additional 
saturation compared to the simpler potentials. This is 
due to 1) a less attractive wjg because much of the bind- 
ing of the deuteron is now provided through tensor cou- 
pling to the channel, and 2) the ability of the tensor 
interaction with third particles to change an attractive 
pair into a repulsive '^P pair This saturation is suffi- 
ciently strong to underbind all the nuclei with respect to 
experiment, and it opens the A — 8 mass gap by making 
^Be less than twice as bound as ^He. However, it leaves 
the A = 6,7 nuclei with only marginal stability. The ten- 
sor forces begin to mix the ^■^+^i[n] states so they are 
no longer eigenstates, but several sets of states, like the 
2+ and 3^ states in ^Li, remain nearly degenerate. 



The spin-orbit terms in AV8' provide much more mix- 
ing and clearly break the J'^ degeneracy, producing a 
spectrum that is properly ordered in the ^ < 8 nuclei, 
although the splittings of most spin-orbit partners are 
smaller than observed. The binding energies shift slightly 
compared to AV6', some up and some down, with the 
A = 6,7 nuclei becoming more stable while the A = 5,8 
mass gaps are preserved. Going to the full AVIS inter- 
action decreases the overall binding energies slightly, be- 
cause the less attractive higher partial waves in AW scat- 
tering are now accurately represented, but the relative 
excitations are virtually unchanged. However, the energy 
differences among isobaric multiplets, like the ^Li,^Be 
mirror nuclei, are improved due to the addition of the 
extra EM, CD, and CSB terms in AV18. 

Finally, the addition of the IL2 NNN potential to AV18 
gives a very accurate description of the light nuclei. It 
adds the necessary additional binding that the realistic 
AW potentials lack and increases the splittings among 
spin-orbit partners. It makes ^He and '^He strong-stable, 
and it produces the correct 3+ ground state for ^''B, 
where all the simpler models (and other realistic WW in- 
teractions like CD Bonn |^^) incorrectly give a 1+ state. 

Many of the results from these models can be under- 
stood by simply counting the number of ST pairs, Nst, 
in the IP wave function for a given state and applying 
a weight factor appropriate for the potential. A good 
estimate of the relative binding for the AV4' and higher 
models is obtained using weights from the OPE poten- 
tial: 3Wio -t-3Woi — Wii — 9Woo. This reflects the approxi- 



4 




FIG. 2: Nuclear energy levels for the more realistic potential models; shading denotes Monte Carlo statistical errors. 



mately equal attraction in '^S and potentials, the small 
repulsion in "^P, and the large repulsion in ^P. The values 
for ^He, ^He, ^He, ^Li, ^Li, ^He, ^Be, and i°B ground 
states are 18, 18, 21, 21, 27, 24, 36, and 39, respectively. 
With this estimate, the marginal stability of ^He and *Be 
against breakup and the roughly equal binding of ^He 
and ^Li are expected. It also provides the order of ex- 
cited states and thus the relative amount of mixing, e.g., 
in ■^Li the 2p[3], 4p[21], and 2p[21] states get weights of 
27, 21, and 15, as the number of S-wave pairs decreases 
going from [3] to [21] symmetry, and the ratio of '^P to 
^P pairs decreases going from quartet to doublet spin. 

For AV6' and up, the importance of the OPE potential 
is evident from its expectation value, which is typically 
80% of (vij) 1^. These findings are consistent with the 
important role of the spin-isospin interaction in fixing the 
shell gaps in nuclei ||l|, and support a close connection 
between nuclear structure and the underlying features 
of QCD, particularly the special role of the pion as the 
Goldstone boson, and the dominance of spin-isospin and 
tensor forces in 1/Nc expansions |^]. 

We see from the present studies that purely central 
nuclear forces are nonsense for nuclei beyond the s-shell, 
where it is crucial to incorporate the difference between 
attractive even and repulsive odd partial waves. While a 
mo del like AV4^ 



V. R. Pandharipande for many useful discussions. The 
many-body calculations were performed on the parallel 
computers of the Mathematics and Computer Science Di- 
vision, Argonne National Laboratory and the National 
Energy Research Scientific Computing Center. This work 
is supported by the U. S. Department of Energy, Nuclear 
Physics Division, under contract No. W-31-109-ENG-38. 



can produce the eneigy satmatiuu and 
clus tering that appears in the p-shell, our model calcu- 
lations suggest that obtaining the mass gaps at ^ = 5, 8 
and stable A = 6, 7 nuclei is a very sensitive issue, and 
may well require both tensor and spin-orbit forces as in 
the AV8' model. Finally, to get a truly good fit both to 
the ground state binding energies, the spin-orbit split- 
tings in the excitation spectra, and (in the case of ^"B) 
the ordering of spin states, we need multinucleon forces. 
We thank J. L. Friar, D. Kurath, D. J. Millener, and 
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